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1 Introduction 



Interesting characteristic features of a given shape variable can be appreciated only in calcu- 
lations to next-to-leading order. One of the most interesting variables is the jet-broadening 
^-distribution, introduced in Ref. [1], in which 2B is the sum of the moduli of the transverse 
momenta of all emitted particles with respect to the thrust axis in units of Q, the total final 
state mass. 

Consider, in perturbative QCD, the radiation emitted in e + e~ annihilation consisting of the 
primary quark and antiquark with 4- momenta p and p and the secondary partons fcj. Hereafter 
we shall consider small values of B so that the secondary partons are soft and the primary p 
and p belong to opposite hemispheres. To be specific we define the right hemisphere as the one 
containing the quark p. 

The simplest variable is the single-jet broadening (the right-jet broadening Br in what 
follows) defined by 

2B R = Y)ki\ + Pt, Pt = ~^k ti . (1.1) 

Since the transverse momenta are taken with respect to the thrust axis, the total vector sum 
of transverse momenta in each hemisphere is zero. One introduces the total jet-broadening Bt, 
as the sum of right- and left-jet broadening (Bt = Br + Bl), and the wide-jet broadening 
B w = Max{B R ,B L }. 

We consider first the ^-distribution. The jet-broadening distribution for the right hemi- 
sphere is given in terms of the multi-parton emission distribution da n by 

da d UBr), 



a(d\nB R ) d\nB R 

(1.2) 



Ir(Br) = W ^ e Ubr - Y, \hi\ - 

n J a \ ieR 



ieR 



For small Br, which corresponds to all the final-state partons having small transverse momenta, 
one can approximate da n as a product of two factors. The first is a coefficient factor which 
depends only on a s (Q). The second is an evolutionary exponent, which describes the production 
of small- kt partons off the primary quark-antiquark pair. Hence 

I R (B) = C(a s (Q))X R (B,a s ). (1.3) 

The essential momentum scales in the coupling in T>r range from BrQ to Q. The perturbative 
treatment that we shall pursue requires BrQ 3> Aqcd- The accuracy of the perturbative 
treatment is limited by non-perturbative power corrections of relative order (In Q)Aqcd/(BQ) 
which are treated elsewhere [2, 3]. 

If only one gluon is present then the quark and gluon transverse momenta are equal and 
opposite — hence the natural factor of two in the definition of Br (1.1). But starting from two 
gluon emission, the situation already becomes significantly more complex, because the modulus 
of the quark transverse momentum depends not only on the moduli of the gluon transverse 
momenta but also on their relative angles. So in higher orders one has to control simultaneously 
the sum of the moduli of the momenta, and the modulus of their vector sum. It turns out that 
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in the kinematical region where a s lnl/BR <C 1, one gluon has a transverse momentum much 
larger than that of the others, so that the quark contribution is still ~ 5 -Br and therefore 
easily accounted for [4]. On the other hand, in the region of extremely small B, a s In 1/Br > 1, 
the gluonic contribution to Br comes from several gluons with comparable transverse momenta, 
so that the problem of an accurate treatment of the quark recoil becomes severe. 

A perturbative analysis based on the all-order resummation of leading and next-to-leading 
logarithmically enhanced contributions to jet-broadening distributions was performed by Catani, 
Turnock and Webber (CTW) in [4]. Their result in the kinematical region B <C 1 is based on 
the "exponentiation" of one gluon emission. The approach that they developed guarantees 
that in the corresponding exponent, all terms with a™ ln n 1/B,m < n are kept track of. At a 
given point, they made the approximation that the quark recoil contribution to Br is always 
= \Br. This simplifies the answer, but mistreats the a^\a n l/B terms, starting from 
n = 2. In this note we will show how to treat properly the quark recoil and thus how to improve 
the CTW prediction for the jet-broadening spectrum. 

Before discussing our result we first recall the form of the double-logarithmic contribution. 
Here one can simplify the analysis by assuming \pt = \Br- Moreover one assumes that all 
the final-state partons are soft and collinear gluons, emitted independently. This small Br 
contribution to the jet-broadening distribution is given by 

Zr(Br) = e -flo(VB«) (1 + O {a n ln n 1/Br)) (l4) 

where Rq(1/B) is the soft part of the gluon emission distribution (see Appendix A) integrated 
over the region kt > Br 



m/BR) _ f f r dz « 1 _ ^sel _> 1/Br + . . . (L5) 

Jq 2 b 2 r «t Jk t 2vr vr 



The limitation k t > Br comes from the fact that for small Br, to leading order the real emission 
takes place only for kt < Br and here is cancelled by part of the virtual contribution. Thus only 
the virtual contribution remains in the region kt > Br. Taking into account the running of the 
coupling, i?o becomes a series with terms of the form a™ ln™ +1 B, while the neglected terms are 
one power of In B down. 

In our analysis, we intend to compute all single-logarithmic corrections, i.e. corrections of 
order a™ ln ra B. We show that to achieve this accuracy it suffices to "exponentiate" the next-to- 
leading order single-gluon emission formula [4] and to treat properly the quark recoil. The final 
result can be written in the form: 



T e {1 _ R/) J e - fl (w (1 + 0(<ln n-i 1/SH ) )i n > 2) 



_ dR(l/B R ) 
dlnB R 



(1.6) 



where the radiator R(\/2Br) is given by the one-gluon emission distribution, as in (1.5), but 
obtained from the next-to-leading order splitting function. The argument of the radiator is the 
jet-broadening Br rescaled by a function 2/A which has a single- logarithmic expansion, i.e. with 
leading terms of the order a™ ln n 1/Br. For Br finite, A tends to 2, while for small Br it tends 
to 1. The function A takes into account the effect of quark recoil. If one neglects quark recoil 
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and puts for instance the quark transverse momentum at p t = Br, as done in [4], then one finds 
exactly the same expression (1.6) but with the function A frozen at A = 2. 

Within our calculation we are also able to show the absence of the first non-logarithmic 
correction to 1(B) that comes from the running coupling at low scales, a s (BQ). Though for- 
mally classified as being subleading, O (a"ln n_1 1/-B), numerically it could be dangerous, since 
a s (BQ) increases at small B. 

It may be surprising that the accuracy in (1.6) can be achieved by using an independent 
multi-gluon emission distribution. The important point to realise is that subsequent gluon 
decay can be neglected to single-logarithmic order. This will be discussed in detail in the 
paper. However the following simple argument explains why. The i?-spectrum contains a 
characteristic exponent of a s In 2 1/Br. A value of B can be changed by a non-collinear non-soft 
gluon decay, with the transverse momenta of both offspring partons being O (BQ). Given the 
relative probability of such a decay, O (a s ), we get a correction B — > B(l + a s ), which translates 
into a s ln 2 (B(l + a s )) = a s ln 2 B + a 2 In 5, the latter being a negligible effect with a power of 
the log-enhancement factor smaller than that of a s . 

The paper is organised as follows. In Sect. 2 we derive the Br distribution using the 
independent emission distribution and by including the effect of quark recoil. We generalise 
the analysis to the total and wide jet-broadening distributions Y^t(B) and Ejy(-B). 

In Sect. 3 we show that our result can be derived from the CTW equation obtained from 
the coherent branching process. To achieve the desired single- logarithmic accuracy one has to 
improve the treatment of recoil. We also verify that one can neglect the gluon branching process, 
thus justifying the use of the independent gluon emission approximation. 

In Sect. 4 we perform the actual evaluation of E(-Br) to the required accuracy. 

In Sect. 5 we present the results of numerical analysis. We compare the final result with 
O (a 2 ) numerical results obtained with the EVENT2 program [5]. We give all the necessary 
information to be able to carry out the matching with the fixed order perturbative results. 
We show the difference between our result and that of CTW, and compare the matched and 
non-matched calculations. 

In Sect. 6 we discuss, comment and look forward. 



2 Independent gluon emission 

To compute jet broadening to single-logarithmic accuracy, it is enough to use the independent 
gluon emission distribution from the primary quark-antiquark pair. This is given in terms of 
the q — > q + g splitting function to next-to-leading accuracy. 

The independent gluon emission distribution is 

1 n d 2 k r 

dw n = — Y[ —r j-dzi 2P qq [a s ,Zi] Q(zi - k u ) V , (2.1) 

where the factor 2 takes into account the fact that, for each gluon, the integration domain 
includes both the right and left hemispheres. The quark splitting function is given, in the MS 
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scheme by 



P qq[ a MS> Z } = C F 



a m( k t) 1 + z2 



2vr 



1 + 



a 



MS 
2vr 



K 



+ 



18 6 



(2.2) 
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for rif flavours. The non-soft part of the quark splitting function, 1 — z ~ 1, produces a single- 
logarithmic contribution to the exponent, and must be kept, while the neglected part of the 
two- loop anomalous dimension generates negligible terms of order a 2 InB. The phase space 
region is (1 — z) > kt/Q (see Appendix A). 

We work in the physical (CMW) scheme [6] in which a s is the measure of the intensity of 
soft emission. We then have 



P«[°., z] = ^ , a s (k t ) = a m (k t ) (l + ^ K + O (a 2 m ) ) 



To this order, the virtual correction V is given by 



InV 



^ dk 2 r 1 -^ 



f dK f 
J % Jo 



dz 2P qq [a s ,z] 



(2.3) 



(2.4) 



For collinear-safe inclusive quantities one can integrate over kt in (2.4) and ku in (2.1) down to 
zero. 

The distribution (2.1) is normalised to unity: Yin I d-Wn = 1; where the transverse momen- 
tum of each real gluon is integrated up to Q. The independent emission distribution is valid 
only for small ku- It simplifies the treatment but mistreats the non- logarithmic region of large 
transverse momenta, kt ~ Q, both in real and virtual terms. This is compensated by the factor 
C(a s ) in (1.3). For small B values this factor is S-independent. The independence is embodied 
into the S factor, 



Efl(B) = Y, [ dw n e(2B-Y^ \h 

n J \ ieR 



ieR 



(2.5) 



To make use of the factorisation structure of the multi-gluon matrix element (2.1) we introduce 
the Mellin representation of the ©-function 



el2B-J2\ku\ 

V ieR 



J2hi 

ieR 



du_ e2uB d 2 Pt d 2 b e - ift .b e -„ pt 



n 

ieR 



e -vk t i e -ib-kti 



(2.6) 



Here v is the Mellin variable conjugate to B, which runs parallel to the imaginary axis, to the 
right of v = 0. We have introduced the integration over the quark transverse momentum p t and 
the constraint Pt = — Y2ieR ^ti is implemented by the integration over b. We obtain 



dv 
2ttiv' 



-,2uB 



a{v) 



where the contour lies to the right of all singularities of a(u) which is given by 



d 2 Pt d 2 b 



e -m-b e -vpt . e ~n(u,b) 



vbdb 



o [y 2 + b 2 fl 2 



-K(v,b) 



(2.7) 



(2.8) 
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with the radiator 



H(v,b) = [ ^dzP[a s ,z] \l - e - ukt e- iUt ] . (2.9) 
J Trkf L J 

In the CMW scheme the two-loop a|? contribution to the quark splitting function (2.3) is regular 
at z = 1 and can be neglected because it produces corrections to (2.5) of order a™ In™ -1 1/B. In 
the Ms regularisation scheme instead one should explicitly keep the infrared-singular contribution 
of order to the splitting function (2.2), in order to reach the desired accuracy, as done in [4]. 

Notice that one is completely inclusive with respect to gluons emitted in the left hemisphere. 
Thus their contribution cancels y/V, i.e. half of the virtual corrections. The remaining 
corresponds to the term 1 in the square bracket of (2.9). 

Performing the azimuthal integration we obtain 

K(u, b)= 4(1 /h) [l - e-" fe Mbh)} , (2.10) 

Jo fc t L J 

with (j>(l/kt) the one-gluon radiation formula to next-to-leading accuracy 

cf>(l/k t ) = 2^M^ (in I - |) + O («*(*,)) . (2.11) 

In the leading double-logarithmic approximation Tir(B) is obtained as follows. Since quark 
recoil is irrelevant to this order, we are free to set pt = -Br, which corresponds to neglecting the 
6-dependence in the radiator, i.e. replacing lZ(v,b) with lZ(v, 0) in (2.8). Then, for large u, one 
approximates [1 — e~ ukt ] — > &(k t — \/v) and further approximates the ^-integration by setting 
v — > 1/B. This leads to the leading order result (1.4). 

Following CTW, from a(v) one deduces the total jet broadening Er(-B) and the wide-jet 
broadening E^/(-B) distributions given, for small B, by 

Et(*) = / ^-e 2 - , = ( [ ^-e^ a^Y . (2.12) 



2-Kiv " ' \J 2mv 

Characteristic values of the Laplace parameter v and of the impact parameter b that determine 
the final answer for the ^-distribution satisfy v ■ BQ ~ 1 , b ~ u, hence, in the B <C 1 
kinematical region we have a large parameter vQ ~ 1/B 3> 1. 

A systematic evaluation of both cr{y) and S(i?) to the required accuracy will be considered 
later. First we discuss the connection to the approach of CTW, obtained in the framework of 
coherent branching. 



3 Coherent branching 

In [4] a technique for analysing T,[B) was developed, based on the evolution equation for the 
distribution T q (Q, k, P t ) of the variable 



1^1 + 


^2 hi 
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in a quark jet produced with vector transverse momentum k. The final physical distribution is 
obtained by setting k = 0. The distribution a(v) in (2.8) is given by its Laplace transform at 
k = 

poo 

a(v) = f q (Q,6,v) , f q (Q,k,u) = dP t e^ k - p ^T q (Q,k,P t ). (3.1) 

Jk 

From the coherent branching picture one has the following evolution equation (see [4]) 

f q (Q, ^v) = l + j Q Q ^J Q 1 dz P qq [a s , z] {^(|S|-l**l-ld-.)S-«l) . 

f q (zq, zk + q t ; v)f g ({l - z)q, (1 - z)k - qf, v) - f q (q, k; !/)} , (3.2) 

where qt = z(l — z)q, qt = \qt\ and k = \k\. Here, T g is the corresponding distribution for a 
gluon jet. The first term in the curly brackets describes real parton splitting, while the second 
subtraction term accounts for the virtual effects. 
Introducing the Fourier transform 

T i (Q,b,u) = J^e^- vk T i (Q,k^), i = Q,9, (3-3) 

we find 

r Q2 da 2 r 1 

T q (Q,b,v) = N(v,b)+ 4- dzP qq [a s ,z] 

Jo Q Jo 

{/ ^Mrz(l- z)q)T q (zq,b+(l- z)r,u)T g ({l-z)q,b- zf,u) - r,(g,6,i/)|, (3.4) 
where N(u, b) is the Fourier transform of the inhomogeneous term 

Since the dependence on the evolution parameter Q is contained in the upper limit of the 
g-integration, we consider the logarithmic derivative of (3.4): 



r;(Q,M = ^r,(Q,^). (3.6) 



Then q in the integrand gets replaced by Q and q t = z(l — z)Q. It is straightforward to verify 
that the integral term possesses the damping factors e ^ uz ( 1 - z )Q J (bz(l — z)Q) , and therefore 
is concentrated at small values of (1 — z), 

(l- Z )<mm{(uQ)-\(bQ)- 1 }. (3.7) 

As a result, any correction proportional to (1 — z), or Inz, will produce a power-suppressed 
contribution ~ 1/uQ oc B <C 1. Taking into account that r ~ b ~ u, this allows us to 
approximate q t = z(l — z)Q ~ (1 — z)Q, which correction is of the order of a s (Q)/7r and will 
be absorbed into the coefficient function C(a s ) in (1.3). We can also replace in (3.4) 

T q (zQ,b+(l-z)r,v) => T q {Q,b,v), 
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and 

T g ((l-z)Q,b-zr,u) => T g {(l-z)Q,b-r,v). 

The accuracy of these approximations is of the order of the relative corrections a s /ir • B. Notice 
that the first argument in r g , generally speaking, cannot be expanded in (1 — z) <C 1 since the 
corresponding dependence is double-logarithmic. 

By making these simplifications (qt = (1 — z)Q) we obtain 



r,(Q,6, 

The final solution reads 

72 



^ = 2 / >1 ^P ?? [a s ,z]{Jo(6 9t )e-^r ff ( % , % ,^)-l}. (3.8) 
V) Jo 



f q (Q,0,v) = J ^N{u,b)e 



(3.9) 



S(u,b)= -| / dzP w M<ft)>z] 1 - e-^ + ^T 9 (g t , % ,^) 

JO JO L 

Neglecting gluon branching, T g ==> 1, we obtain the result of previous section 

<7(i/) = f,(Q,0;i/) = (2.8). (3.10) 

To see that gluon branching, T g ^ 1, is indeed negligible within our accuracy, we look at the 
correction to the exponent in (3.9): 

55(^6) = / -4 / dzi^Mft),*] e~^+^ [T g (q t ,q t ,u) - 1] . (3.11) 

Jo Kq t J 

First we notice that the exponential factor forces qt < v~ x . On the other hand, by examining 
the evolution equation (3.2) we observe that T g (q t , q t , u) — 1 vanishes for q t <C v~ x . As a result, 
the (^-integral is concentrated at q t v ~ 1, T g (q t ,q t ,v) — 1 ~ q s £(ln(<ftz/)) , and the correction 
amounts to 

5S ~ a s 2 ln(K?) , 

with the single logarithmic factor emerging from the z-integration. 

This first perturbative correction is due to quark — > (quark + two gluons/gc/) pair splitting 
processes in which two secondary partons have similar emission angles and energies of the same 
order. It is easy to check that a correction of the same order originates from non-collinear 
two-parton production at large angles ~ 1, and, in particular, of the configuration of partons 
falling into opposite hemispheres. In the present treatment the left and right jets contribute 
independently to the event broadening and so such subleading contributions are not included. 



4 Evaluation of E(B) 

In this section we first evaluate the radiator TZ(v, b) and then perform the 6-integral in (2.8) to 
evaluate cr(v). For the sake of simplicity in this section we put Q = 1. 
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4.1 Radiator 

To evaluate the radiator for large v and b values we introduce the function 



R(p) 



/ 1 rlk / 
^(k- 1 ), » = h(v + V^+V), (4-1) 



which corresponds to (1.5) with the next-to-leading splitting function. With account of the 
running coupling, it contains all terms of order a™ ln n+1 /i as well some of the essential a™ ln n /i 
terms. We write 

K(u,b) = R(p) + 5R(v,b) , 

and study the correction 

rlk / \ f 1 rlk 

SR{y, b) = ^ ^{k' 1 ) (l - e' kv J (bk)) ~ J i "f ^ k ~ l ) e~ kv J {bk) . (4.2) 

It is determined by the integration region k ~ /i -1 and can be written as (e — ► 0) 

f 1 rl7 f'°° rl7 

SR(v,b) = / —z e rb(p/z)- —z e e- zu ^Jo(bz/fi)rb(fi/z) 
Jo z Jo z 

= -Sfa/z)- -/e^^)- -e-'^iMbz/ri-WWz), (4.3) 

JO z JO z JO z 

where we have neglected contributions of order of exp(— fi). 
To evaluate 5R we expand (f)(fi/z) in powers of In z 

rb{ii/z) = BUji) - R"{n) In z + . . . 

with 

aln/x 

= -^(m) = 2C,^iM + O (a? ln M ) , ^ 
aln/x 7r y 

R'"(n) = O (a s 2 ln/i) . 

We have that R ! "{[i) is beyond our accuracy since it would lead to corrections to T,(B) of order 
a™ ln n_1 1/B. In R"(n) we keep the correction of the order of the coupling at the reduced scale 
a s (l/fi) and neglect its higher powers in a s . To our accuracy we could neglect the scale 1/fi 
in a s . However we shall keep track of the non-logarithmic corrections proportional to a s (l//i) 
that emerge in the course of approximate evaluation of the b- and z^-integrals. This enables us 
to guarantee that no first-order correction a s (BQ) is present in the final answer. 

The contribution to 5R(v, b) proportional to R'(n), the leading order correction, is given by 

5RW = #( M ) • lE . (4.5) 

Having chosen a different definition of fi would have resulted in an additional logarithmic con- 
tribution. 

The non-logarithmic correction, of order R"(n) ~ a s (l//x), is 

5RW = R "^) . A(v, b) , (4.6) 
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(4.7) 



where (see Appendix B) 

A(v,b) = U (in ^ + 7E ) 2 + \ + In /i In ^ + c(i/, 

c(i/,6)= / — Inx e~ ux ( J (bx) - 1) , 
Jo x 

We have then 

H(u, b) = R(n) + SR(u, b) = R{^ E ) + R"{n) • (A(i/, 6) - ± 7 I) + O (#"') . (4.8) 
Expanding in In (i/v we obtain 

b) = R(p) + In ^ + R"(u) ■ A + O (#"') . (4.9) 

where 

P = z/e 7B , (4.10) 

and 

A(/x/i/) = ±V'(1) + In £ - In v In £ + c(v, b) , (4.11) 
is a function of the ratio n/v and has no lnv terms (see Appendix B). 

4.2 a(u) 

Now we substitute (4.9) into (2.8) and performing the b-integral to evaluate a(v). Changing the 
integration variable to y = \/v 2 + b 2 we obtain (/x/V=|(l + y)) 

*>dy (l + yY R ' {P) / , + y 



To calculate the main contribution we introduce the rescaling function A(i?') 

/,.- P '/// /i • u s 



w ""-y, #w • (413) 

To estimate the correction R" ~ a s to first order we can neglect the exponent R' since A is of 
order one (it has no Inu contribution). It equals (see Appendix B) 

/°°^A(i±^)=I(^(l) + ln 2 2). (4.14) 

Finally we arrive at 

a(u) = e~ R W ( A )-*M [i _ B»(p) i (V /(1) + In 2 2) + O (R'")] 

= e ~R(uX) ^ _ i R »^ _ i jR "(p)(i n 2 2 - in 2 A) + O (R" f )] , ' ) 

Since A is a function of R'(P), the leading terms of its expansion are a™ In™ v. For small values 
of R'(v), the function X(R') is close to 2, 

X(R') = 2- R'{u) " 2 In 2 2^ + O (i? /2 ) , (4.16) 
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while for large values of R'(v) it approaches unity: 

Perturbatively, A is close to 2, see (4.16). Therefore the quantity R"(ln 2 X — In 2 2) in (4.15) is 
of order R'" and thus negligible. We can write 

a (y) = e ~ R ^ [1 - 1/(1) R"(v) + O (R'")] . (4.18) 

Notice that this estimate is uniform in R'. Indeed, for large values of R' the In 2 A term in (4.15) 
vanishes, see (4.17). At the same time, the In 2 2 term disappears as well, since for large R' the 
integral (4.12) for the correction is concentrated near y = 1 where A(l) = instead of 

(4.14). 

The next-to-next-leading correction R"(p) ~ a s {Q/v) ~ a s (BQ) that we kept in (4.18) will 
cancel, in the first order, against a similar correction coming from evaluation of the ^-integral. 

4.3 The ^-integral 

Here we finally compute the integral over v in (2.7) to obtain To this end we expand 

o(y) around some point vq in powers of \o.{v/uq) = Int. Introducing rjo = i?o\(R') we have, to 
the required accuracy 

Efl(B) = e- 11 ^) f — e 2 » QBt t^'^ [l - \R"(u) (in 2 t + y/(l)) + O {R'")] . (4.19) 

J 27T7 t 

Notice that the corrections coming from differentiation of X(R') in the argument of R(v\) do 
not contribute: A ~ (a s lni>) n , therefore R' ■ A' ~ a™ln m_1 u, escaping our resolution. 

Now we have to choose the value of vq such as to keep the characteristic value of lni not 
large. The leading term is given by the following basic integral 

27Tit e t r(l + 7 )' 7 {V0h [ ' 

and we obtain 



-R{vo) 



! ~ + + ° H ff^y ' (421) 

The correction proportional to R"{r]o) becomes 

(3? + *' (1> ) irS = O 1 "* 2 ^ 1 - * (1 + ]2 + - + ww) ■ (4 - 22) 

Choice of u . To ensure the smallness of the correction one has to choose vq in such a way 
that o\{2Buq) remains finite. By choosing 

= , = H>e» = ^ , (4.23) 

the first correction of order R" in (4.19), given by (4.22), becomes uniformly small: 



1 nil 

2 K 



Ml)-^(l+ 7 )) 2 + V/(l)-^(l+ 7 ) • (4.24) 
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It vanishes in the first order, that is for 7 = 0. 
The parameter rj reads 



and we arrive at 



Vo = ^0 A = ^ B 1 , 

£*(B) = e^ (A/2B) j^q^y - 4 0) ( 25 /A) ■ (4-25) 

This first form of the answer is the CTW distribution with the important modification that the 
argument of the double- logarithmic radiator is 2B / A instead of simply B. Since R' enters only 
in the prefactor, we can take R' = R'(l/B). 

We can further simplify the answer by expanding R(\/2B) in powers of In A/2. Neglecting 
the contribution O (R"), we arrive at the second form of the answer, 



where R' = R'(l/B). This is the CTW answer modified by the single-logarithmic factor (2/\) R> . 
The two forms of our answer, (4.25) and (4.26) are both correct (and equivalent) to single- 
logarithmic accuracy. 

Finally our results for the total and wide jet broadenings (given in a form analogous to 
(4.26)) are 

/ 9 n 2R> p -2 7B i?' / 9 \ 2R' 

^-(! -WtW) • ^ (B) , (4,7) 



and 



9 x 2R' p -2 lE R! / o \ 2R> 



respectively. As before, i?' = R'(l/B). The value of C(a s ) in (1.3) is 

C(q s ) = 1 + Y, C n a£ , C ^ = % [f ~ y) • ( 4 - 29 ) 

n=l ^ ' 

For consistency with the order at which we have performed the resummation, it is sufficient to 
know C only to first order in a s . 

5 Numerical analysis 

5.1 Comparison with two-loop result 

We compare the single-logarithmic result for the total and wide jet distribution [i = T, W) 

doj _ dI t (B) 
a(d\nB) dlnB ' 

with the exact two- loop distribution obtained with the EVENT2 program [5]. 

To obtain the contribution from our calculation one needs to expand eqs. (4.27,4.28) up 
to second order in a s and multiply by (4.29). At order a? the result should accurately reproduce 
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the coefficient of terms a 2 ln m 1/B with m = 2,3, 4. To verify this, in fig. 1, we plot as a function 
of B the differences 

, s -iog rm = ^Bg) dir lo9 (B) 

1 K >~ dlnB dlnB ' 



s *-io g(B) = dlf(B) _ dlt^(B) 
i K ' dlnB dlnB ' 

with lf{B) obtained from EVENT2, lf~ l ° s (B) obtained from the expansion of the resummed 
result to double-logarithmic accuracy (a 2 ln m 1/B with m = 3,4) and I- g (B) to single loga- 
rithmic accuracy (a 2 ln m 1/B with m = 2, 3, 4). The comparison is made within the MS scheme, 
and with / normalised to the Born cross section rather than the total cross section (as this is 
what is supplied by the EVENT2 program). We show only the coefficient of the (a s C^/27r) 2 
part of 5i, as this is the only component modified by our new treatment of quark recoil. 
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Figure 1: The coefficient of the (a s CV /2n) 2 component of the difference, Si(Bi), between the numerical 
two-loop calculation of the ^-distribution, performed with EVENT2, and the expansion to single and 
double logarithmic accuracies of the resummed result. Shown for the total and wide broadenings. 



The quantity Sf l ° 9 (B) is of order a? Inl/B, corresponding to the absence in lf~ log (B) of 
the a 2 In 2 1/B term. Meanwhile, 8 s i ~ l ° 9 (B) is of order a 2 , indicating that in / we have correctly 
taken into account the a 2 In 2 1/B term. 



5.2 Matching 

When comparing with experimental data one often chooses to match the resummed calculation 
with a full second order calculation (as for example from the EVENT2 program). Here we give 
only a brief summary of two main matching procedures (log-R and R matching schemes) and 
refer the reader to [7] for detailed information. 
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For convenience we write 



00 n+i 

n=l m=l 

= Lg 1 (a s L) + g 2 (a s L) + a s g 3 (a s L) + 



(5.1) 
(5.2) 



with L = lial/B. The resummation procedure provides g±, <?2 and C\. Given the resummed 
calculation, the full two- loop calculation contains information on G21 (the first term of g%) 
and C2 and on a "remainder" which does not necessarily exponentiate. Matching schemes put 
together these different parts with the principle ambiguity being the treatment of the remainder: 
the log-i? scheme makes the approximation that it exponentiates, the i?-scheme that it doesn't. 

Here we provide the information that is needed for implementing the matching with our new 
calculation. We recommend the use of the second of the forms in (4.27,4.28), using as the 
equations (18-22) of [4]. The factor by which it is multiplied is 



2R' 



00 dx 
x 2 



1 + x 



-R r 



and one should use 



R' 



2a s C F 



\n(l/B) 



vr l-2a s /3oln(l/B) : 



(5.3) 



(5.4) 



with (3q = (IICa — 2nj)/127r. This particular form ensures that the resummed calculation 
contains only g\ and 52 and no spurious gs type terms, as is required for correct matching. 

In addition the matching schemes require Gi m , G2 m , C\ and C2 explicitly. Gi m , G23 and 
C\ are correctly reproduced in [8] 1 . The new analytic forms for G22 are 



G22 — 



1 



(2vr) 2 

for the total broadening, and 

G22 = 



32 In 2 2 + - vr 2 



C\ + 



2 2 35 

r —9 



C A C F + -C F n f 



(27T) 



-(32 In 2 2)0$ + -w 



35 



C A C F + -C F n f 



(5.5) 



(5.6) 



for the wide broadening. 

Finally for the i?-matching procedure one requires explicit values for G21 and C2. These 
can be determined from the calculation performed with EVENT2. First one determines G21 by 
fitting the a 2 ln(B) component of 1(B); then one subtracts it out and fits the asymptotically 
constant part of 1(B), which yields C2. One obtains: 





Bx 


B w 


G21 


1.988 ±0.25 


1.869 ±0.25 


c 2 


2.330 ± 0.25 


2.946 ± 0.25 



It is necessary to go to extremely small B 



0-8 



before a sufficiently asymptotic behaviour sets 



in. As a result the value of C2 which one obtains is strongly dependent on the value taken for 



Note that the coefficients in [8] are given for (as/27r) n ; this differs from the convention used here. 
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G2i- The results shown for Ci were obtained by fixing G21 (at its central value), and the error 
is that purely from the fitting procedure for C2. Roughly, the value of C2 that would have been 
obtained with a different G21 would be shifted by (—8 ± 2)[G2i — G21 (central)]. 

A final aspect of matching that needs to be mentioned is that one usually introduces a pa- 
rameter By im at which the resummed calculation is constrained to vanish through a substitution 
of the form 

B B B\\ m ^ ^ 

For the purposes of fitting experimental results, a reasonable choice of By mi is important; in the 
section that follows though, it well be set equal to 1 in order to facilitate the comparison with 
results that do not include matching. 



5.3 The B T distribution 

To illustrate the changes introduced by our new calculation, we show in figure 2a the Bt 
distribution in the approximation that A = 2 and in the case with the full treatment of the 
quark recoil. Examining first the case without matching, one sees that the distribution is 
shifted towards lower values of B, as expected. Associated with the shift is an increase in the 
height of the peak of the distribution at small B. 

If one considers instead the results with log-i? matching, one finds that the effect of going 
from A = 2 to full recoil is much smaller 2 . The explanation is linked to the following property 
of log-i? matching procedure: if one supplies it with the wrong forms for G22 and 52, as long 
as those forms are consistent with each other, the input from the full 2-loop result reestablishes 
the correct value for G^- 

In figure 2a one notes that the matched and non-matched curves are quite significantly 
different. At large B this is to be expected, since one is beyond the range of validity of the 
resummed calculation. At small B the difference is less welcome. However there remain certain 
logarithmic exponentiated terms which have not been taken into account in our resummed 
calculation, in particular those that contribute to G21. In figure 2b we plot our non-matched 
resummed calculation with the additional inclusion of the effect of G21 as determined numerically 
in the previous section. At small B this leads to very good agreement with the log-i? matched 
results. At large B, as one expects, a significant difference remains. Figure 2b also shows a 
curve resulting from i?-matching, and it is seen to be everywhere in good agreement with the 
log-i? matching. 



6 Conclusions 

We have shown that a proper treatment of the contribution to B of the primary quarks recoiling 
against an ensemble of soft gluons is essential for predicting the ^-distribution with single- 
logarithmic accuracy. 

2 In contrast it is not consistent to show results for A = 2 with ii-matching, because the ii-matching procedure 
requires the input of values for G21 and C2 — these cannot be reliably obtained if one uses a wrong value of G22- 
The log-_R matching procedure doesn't suffer from this problem because it doesn't require the explicit input of 
G21 and G2. 
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Figure 2: (a) The effect of going from the A = 2 approximation for the quark recoil, to the full treatment 
of quark recoil, both in the case without any matching and in the case of log-i? matching; (b) The effect 
of different matching procedures and a comparison with the case without any matching, but with the 
inclusion of the G21 term; for all curves, y/s = 91.2 GeV and a s (Mz) = 0.118. 

We have verified that the evolution equation for the ^-distribution derived by CTW embod- 
ies all the necessary ingredients to provide the resummed perturbative prediction with single- 
logarithmic accuracy. The improvement we made concerns the solution of this equation. We 
suggested two forms of the final result, (4.25) and (4.26). The former is the CTW-spectrum eval- 
uated at a rescaled value of B — > 2B/X, with A a single-logarithmic function, A = A(a s In 1/B) 
which decreases with B from 2 to 1. The latter form is the CTW-answer supplied with the 
single- logarithmic factor which does the same job of shifting the distribution to smaller B val- 
ues. This tendency is opposite to that expected from the 1/Q power correction effect in the 
.B-spectrum [3]. 

We were able to show the absence of the non-logarithmic contribution proportional to the 
coupling at the reduced momentum scale, a s (BQ), which, if present, could damage the pertur- 
bative prediction at small values of B. 

It should be noticed, however, that beyond the first order in a s (BQ), given the present state 
of the art, such damage looks unavoidable. Indeed, consider the most interesting feature of the 
^-distribution which is its characteristic maximum at B = B max <1. A maximum emerges as 
a result of an interplay between the first-order peaked cross section, oc (a B In 1/B)/B, and the 
all-order Sudakov suppression exponent, exp(— R(X/2B)) ~ exp(— a s In 2 1/B). The latter factor 
takes over, clearly, when R'(l/B) ~ a s (BQ) In 1/B approaches unity. With B decreasing In 1/B 
increases and so does the running coupling a s (BQ). Formally speaking, in perturbation theory, 
that is for Q — ► 00, the expansion parameter a s (BQ) stays small in the region of the maximum. 
However, in reality (and for any foreseeable energies) a s (B max Q) becomes numerically large. 
This undermines the reliability of the perturbative prediction for B < B max since the neglected 
subleading corrections of the order a 2 In 1/5, and among those, a s (BQ)R' ~ a s (BQ) are no 
longer numerically negligible. Corrections of this sort arise, in particular, from the hard-emission 
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subtraction term R' ~ a s In 1/B — > a s (ln 1/B — 3/4) in the single-logarithmic pre-exponent. 

To trigger genuine confinement effects, the (Ai In Q + A 2 )/Q shift in the U-distribution [2, 3], 
it is tempting to look at smaller Q values. However, some care should be exercised here since 
the kinematical range of B shrinks. Moreover, one should bear in mind the above-mentioned 
intrinsic uncertainty of the perturbative prediction for B < B rnax which becomes larger for 
smaller values of Q. To be on a safe side, one should try to stay with the S-values to the right 
of the maximum, B > B max . 

It remains to be seen whether the 1/Q power correction extracted with the use of the 
improved perturbative expression derived in the present paper, will exhibit the expected InQ 
enhancement. 
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A One-gluon emission with single-logarithmic accuracy. 

The single gluon emission distribution is 

do x , , a s C F x\ + x\ 2pQ 2pQ 

=dx 1 dx 2 — — j- r-rj r x i = -^2"> x ^ = ~Fa- 

o 2ir (1 - Xi)(l - x 2 ) Q 2 Q 2 

We introduce the gluon cm. energy fraction z and transverse momentum fc t 

k? (1 - xi)(l - Xl) n 

Q 2 1 - z 

If the gluon is in the right hemisphere we have 

Xl = 1 - \ (z + ^A 2 - 4(1 - z)e 2 ) , x 2 = l-\(z- y/z 2 - 4(1 - z)e 2 ) . 

If the gluon is in the left hemisphere these expressions are interchanged. 
In terms of the gluon variables we have 

^1 = 2 ^dz as ° F 1 + i 1 ~ z ? ~ Zi 1 ~ z ) e2 
(J ~ vr/c 2 Z 2vr ^2 _ 4(1 _ z ) e 2 

where the factor 2 takes into account that the gluon can be emitted in the right or left hemi- 
spheres. Integrating, 



1 , 1 + 1 - z 2 - 2(1 - z e 2 1 3 , 2l ^ /4N 

dz v - 1 v ; = 2 In 4e 2 In e + O (e 4 ) 

y/z 2 - 4(1 - z)e 2 e 2 v ' 



where zq is the zero of the square root and is given by zq = 2e + O (e 2 )- Then for small e we 
can approximate the quark splitting function by replacing 

1 + (1-Z) 2 -2(1-Z)6 2 1 + (1-^)2 h 

V^ 2 - 4(1 - z)e 2 z K Q h 

where the leading correction is of order e 2 lne. It produces a contribution O (a s (Q)) to the 
integral of the gluon emission and is taken care of by the factor C(a s (Q)). 
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B The integral (4.15). 

The coefficient A(i/, b) of i?"(/i) in the expansion of 8R{v, b) in (4.15) is given by 

* V 'J " (B.l) 

1 A A* \ 2 7T 

= A(hi-+ 7£ j + — + In // In - + c(i/, 6) , 

with c(v, b) given by 



z(v,b)= I — lnxe~ ra (J (6x) - 1) =lnz>ln^+ / — lnx e _x (J (xb/u) - 1) . (B.2) 
Jo x v Jo x 



This shows that A(z/, 6) depends only on the ratio /i/V. 
For the relevant integrals over b we use 

vbdb f 00 dy 



^ -ln^= r^lnl±^ = ln2, (B.4) 



(y 2 + & 2 ) 3 / 2 v J 1 y 2 
> vbdb , 2 /j rdy, 2 l + y it 2 



r^kr 2 ^-=r- 2 ^^ = T- in22 - ^ 

7 (i/ 2 + 6 2 ) 3 / 2 f ii f 2 6 
The inverse Fourier integral of c(z/, 6) is 

J {u 2 + h 2 f /2 c ^> h ) = j bdb J Pdpe- up Mbp) ^ — In* e-" x (J (bx) - 1) 

r °° dr 

= — lnx [e~ 2ra -e~ ra ] = - 7 £ln2 + ±ln 2 2 + ln21nz/. (B.6) 
Jo x 

Assembling all terms, we obtain the result in the text. 
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